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ABSTRACT 

We present first results of astrophysically relevant experiments where highly 
supersonic plasma jets are generated via conically convergent flows. The convergent 
flows are created by electrodynamic acceleration of plasma in a conical array of fine 
metallic wires (a modification of the wire array Z-pinch). Stagnation of plasma flow on 
the axis of symmetry forms a standing conical shock effectively collimating the flow in 
the axial direction. This scenario is essentially similar to that discussed by Canto ' and 
collaborators as a purely hydrodynamic mechanism for jet formation in astrophysical 
systems. Experiments using different materials (Al, Fe and W) show that a highly 
supersonic (M ~ 20), well-collimated jet is generated when the radiative cooling rate 
of the plasma is significant. We discuss scaling issues for the experiments and their 
potential use for numerical code verification. The experiments also may allow direct 
exploration of astrophysically relevant issues such as collimation, stability and jet-cloud 
interactions. 
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1. Introduction 



Highly collimated supersonic plasma jets are ubiquitous, occurring in astrophysical 
environments as diverse as Active Galactic Nuclei (AGN, [Leahy 1991 ), Young Stellar Objects 
(YSOs, iReipurth 19"97|) and Planetary Nebulae (PNe, |Soker k Livio 1994| , |Sahai 20001) . In many 
situations these jets represent high mach number flows with remarkable collimation characteristics. 
In the case of Herbig-Haro jets, observed to emanate from newly formed stars, typical velocities 
(V ~ 200 km s~^) and temperatures (T ~ 5000 K) yield mach numbers in excess of M = 20. 
The length to width ratios of the HH jets can be a factor of 100 or more (IReipurth 1997]) . In 
the late stages of stellar evolution (PNe), highly collimated flows are observed which appear to 
change direction with increasing distance from the source ( Borkowski et al. 19971 ). In all cases, the 
ability of astrophysical jets to remain narrow as they propagate over large scales (while changing 
direction) represents just one of many critical questions facing the field. The role of instabilities, 
the effect of jet/ambient-cloud collisions, the collimating effects of magnetic fields and the nature 
of high density knots within the beams are other examples of unresolved issues in jet studies. 

Jets are important not only due to their ubiquity but also represent excellent laboratories 
for the study of important plasma astrophysics processes (shocks, instabilities, etc.). In terms of 
clarifying issues related to the evolution of astronomical objects, the long dynamical or "look-back" 
times, tdyn = Lj/Vj, in jets allow them to be used as proxies for understanding the obscured and 
often unobservable central sources. 

Progress in the study of astrophysical jets is hampered by the complexity of the governing 



(magneto)-gasdynamical equations ( Choudhuri (1998) ). In addition, plasma effects, such as 
ambipolar diffusion ( [Frank et al. 1999 ) or particle acceleration (e.g. Blackman 1996| ) , may become 
important in some instances. Because of the inherent multi-dimensional, time-dependent nature 
of the problem, astrophysical jet studies have relied heavily on the use of sophisticated numerical 



simulations ( Blondin, Fryxell fc Konigl 1990 , 3tone fc Norman 1994 , Cerqueira &: de Gouveia 



Dal Pino 2001 , Gardiner &: Frank (2001) ) . These simulations often include a variety of physical 



processes such as magnetic fields, ionization dynamics, chemistry and radiation transfer. There 
exist few, if any, analytical solutions to the governing equations that can be used as testbeds 
when these processes are coupled together in a multi-dimensional code. In addition, the need to 
simultaneously simulate both large and small scales and the need for fully 3-D evolution means 
many simulation studies remain under-resolved. Thus, in spite of significant progress made by 
simulation studies finding an alternative means for studying high mach number jets, such a 
laboratory experiments, would be highly beneficial to the field. 

Traditional laboratory studies of jets have, unfortunately, been limited to lower mach number 
flows and have not been able to achieve the conditions required to explore key processes such as 
radiative cooling. The advent of high energy devices used for Inertial Confinement Fusion (High 
Power Lasers, Fast Z Pinches) allows for the possibility of direct experimental studies of hypersonic 
flow problems including the affect of microphysical plasma processes. Recent laboratory studies 
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of blast waves, shock-cloud interactions and planetary interiors have shown the promise of these 
kinds of studies ( [Remington et al. 1999| ). High mach number jets have also been the subject of 
these studies ( [Logory et al. 200C , Raga et al. 2001 . Recent work by ( Farely et al. 1999| , jShigemori 



et al. 200q , [Stone et a/."2000D using high power lasers have produced radiative flows whose scalings 



appear to allow contact with astrophysical parameter regimes including radiative jets. 

In this paper we present the first results of a series of experiments designed to study high 
mach number radiative plasma jets using fast z-pinch devices. The purpose of this presentation 
is to introduce and validate the experimental methodology and to provide justification for their 
connection to astrophysically relevant issues. In section 2 we describe the basic experimental 
set-up providing references to the extensive literature on Z-pinch technologies which form the 
basis of our studies. Section 3 provides scaling arguments which show how these experiments can 
be used in addressing astrophysical questions. In section 4 we provide an overview of the results. 
Finally, in section 5 we present our conclusions and discuss the potential for these kinds of studies 
to directly address questions of astrophysical jet formation and propagation. 



2. Experimental Configuration 

The schematic of our experiments is shown in Fig 1. Our studies were conducted using a 
"pulsed power" or "Z-pinch" machine. A review of these types of instruments and their diagnostic 
methods can be found in Ryutov et al. 20"00| and references therein. Our experiments utilize 



plasma flowing off a conical array of wires to produce a directed flow along the array axis. Below 
we briefly introduce and describe the details of the experiment, diagnostics and the component 
physics. 

Acceleration from Wire Array: In our experiments a fast rising current, reaching IMA in 
240 ns, is applied to a conical array of fine metallic wires. The resistive heating of the wires by 
the current rapidly, {5t ~ few nanoseconds), converts the wires into a heterogeneous structure 
with dense, practically neutral cold cores surrounded by a low density hot coronal plasma (see e.g. 
Kalantar and Hammer 1993| ). The global magnetic field generated by the current flowing through 



the array produces the net J x B pinching force which accelerates coronal plasma towards the 
array axis. The characteristic density of the coronal plasma near the wires is n ~ lO^^cm^'^ and 
the measured inward streaming velocity is v ~ 150 km~^ (measured via laser probing diagnostics, 
Lebedev et al. 1999| ). The wire cores act as a reservoir of material, allowing the process of 



formation and sweeping of coronal plasma from the cores to continue for the entire duration of the 
current pulse in the experiment, thus producing a quasi-stationary converging flow of plasma. The 
short rise-time of the current in this system means that current flow is along the path with the 
smallest inductance (at the largest possible radius), remaining predominantly in the vicinity of the 
wire cores. As a result, the streams of the coronal plasma are virtually current-free. Thus in the 
current experiments it appears (based on the value of the current density and MHD simulations) 
that no dynamically significant magnetic fields exist in the coronal plasma flows or the jets they 
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create. 



Conditions on Axis: For purely cylindrical wire array z-pinches, stagnation of the coronal 
plasma flow on the array axis forms a dense, narrow and stable precursor plasma column with 
characteristic density oi p ^ 10^'^ — 10^^ g cm^'^ and electron temperature of T 50 eV. 
The column is confined by the kinetic pressure of the plasma flow (pf ^ ~ pkT, [Lebedev~ 



al. 2001, [Lebedev et al. 1999| ) from the wires. By placing the wires in a conical array, the flow of 
coronal plasma retains an axial component of momentum after convergence on the array axis. The 
conical standing shock that forms on the axis will effectively redirect the flow in the axial direction 
and form a plasma jet. This process has been well studied in astrophysical contexts. Canto and 
collaborators articulated the dynamics of converging conical flows in a series of analytical and 



numerical studies (Canto et al. 1988, Tenorio-Tagle et al. 1988). Those papers demonstrated the 



efficacy of jet production via oblique shocks formed on the axis of the converging flow. In more 



recent works FYank, Balick & Livio 1996 and Mellema & Frank 1997 have demonstrated how these 



converging conical flows can be established in a natural way in the context of stellar wind blown 
bubbles (Young Stellar Objects, Planetary Nebulae). In section 4 we discuss the Canto mechanism 
and its relation to these experiments in more detail. 

Conical Array Experiments and Diagnostics: Our conical wire array experiments were 
performed with 1 cm long arrays composed of 16 fine wires. The small radius of the array was 
8 mm, and wires were inclined at an angle of 30° to the array axis (Fig. 1). Three different wire 
materials were used: Al, stainless steel and W of 25 pm, 25 pm and 18 pm diameter, respectively. 
Parameters of the jets, which were ejected above the array were measured using a laser probing 
system (A = 532 nm, pulse duration 0.4 ns) with an interferometer and two schlieren channels, 
and time-resolved soft x-ray imaging. The interferometer provided two-dimensional data on the 
distribution of the electron density in the jet. The schlieren diagnostic is sensitive to the gradients 
of the refractive index, and the probing laser beam in one of the schlieren channels was delayed 
by 30ns, which allowed us to measure velocity of the jet propagation in the same experiment. A 
Four-frame gated soft X-ray pinhole camera (2 ns gate with 9 ns separation) was used to record 
two-dimensional time-resolved images of the jet emission. Spatial resolution for the laser probing 
was Sx ^ 0.1 mm and for the soft x-ray imaging 6x ~ 0.3 mm. 



3. Scaling Relations 

The purpose of our experiments is to create jets whose properties can be scaled to those 
of real astrophysical systems. Solutions to the hydrodynamic equations which govern jets 
produced in these studies can be shown to be invariant under certain scaling relations ( |Landau| 
fc Lifshitz 1987 ). In general, when dissipative processes can be ignored it is possible to identify 
key dimensionless parameters which will control the evolution of the system. When the scaling 
relations are obtained, two systems with the same values of these dimensionless parameters will 
behave similarly regardless their differences in physical size and temporal duration. The allowed 
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scaling transformations between laboratory experiments and astrophysical systems has been 
articulated by Ryutov et al. 1999] and [Ryutov et al. 2001 



In our jet studies there are 6 dimensionless parameters we need to consider. The first three 
relate to global properties of the jets ( [Blondin, Fryxell fc Konigl 1990D . They are: the sonic mach 
number M = Vj/c; the jet to ambient density ratio r/ = Pj/ Pa] the ratio of cooling length to jet 
radius x = dc/vj. Typical ranges of values obtained in stellar jets are M > 10, r/ > 1 and x ^ 1- 
As we shall show in the description of results, our experiments are comfortably within these 
ranges. 

The second set of dimensionless parameters concerns the microphysical properties of the 
flow. The first is the localization parameter, which is a measure of the fluid- like nature of 
the plasma relative to the scale of the jet radius: 5\\ = Xmfp,\\/fj where \mfp,\\ is the mean 
free path of the particles parallel to the jet radius. The second and third parameters concern 
the role of the dissipative processes: viscosity via the Reynolds number Re = VjVj/v] heat 
conduction via the Peclet number Pe = rjVj/vh (note Pe measures the importance of heat 
convection in comparison with conduction). In these expressions v and Vh are the viscosity 
and thermal diffusivity respectively. Using typical conditions found in the jet beams in 
our experiments (downstream of the conical shock where the plasma streams converge at 
the axis: Vj ~ 200 km/s, rj ~ .5 mm, "/ig ~ 10^^ cm~^,Z '^10, T ~ 50 eV) we find 
S\\ ^ 10~^, Re > W^, Pe > 10. (In determining for the jet, we used the temperature since 
the supersonic bulk flow is perpendicular to the jet radius in the outflow. If instead we consider 
5|| associated with the converging flows toward the axis which initially produces the jet, we must 
then use the bulk speed since that corresponds to the relevant velocity parallel to the jet radius. 
This would then give 5|| ~ 0.1 — 0.3 for this pre-jet converging flow.) Note also that while Pe > 1 
as it should be, the value obtained is low enough such that heat conduction may play a role in 
smoothing small scale features. 

Thus, in general, both the localization and Reynolds number are in the correct regime for the 
jet flow. We conclude that, with caveat that heat conduction effects must be explored in more 
detail, our experiments provide conditions which allow scalings to astrophysically appropriate 
environments, with the caveat that the magnetic flelds are evidently not dynamically important 
inside the jets. 



4. Results 

Fig 2 shows results from laser probing measurements of a jet formed in the tungsten wire 
array. These images show the high degree of flow collimation obtained in case. The jet has a 
sharp, well-deflned boundary and propagates with a velocity signiflcantly higher than its radial 
expansion. The schlieren images of the jet shown in Fig 2 were obtained in the same experiment at 
two different times, and the measured displacement of the jet tip gives velocity Vz ^ 200 km s~^. 
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Comparison of the two images indicates a far smaller velocity of the radial expansion of the jet, 
Vr < 7 km s~^. This allows an estimate the internal Mach number of the jet as M ~ Vz/Vr ~ 30. 
Abel inversion of fringe shifts in the interferogram shown in Fig. 2 yields typical electron densities 
at the jet boundary n ~ 10^^ cm~^. The electron density in the jet decreases towards the axis. 
This, however, does not necessarily mean that the mass density distribution is also hollow. It is 
possible that the decrease of the electron density near the axis reflects the difference in ionization 
levels. Note that since the jets propagate into a vacuum in these experiments we have density 
ratio 77 >> 1 

Variation of the wire material in the experiments allows us to examine the importance of 
radiation losses on jet collimation by converging conical flows. Previous experiments ( [Lebedev ei| 
al. 2000] , Lebedev et al. 1999) with cylindrical arrays indicate that the parameters of the plasma 



flow (an inward velocity and the mass flux) are essentially insensitive to the material used. The 
increase of the atomic number (A), however, leads to a significant increase in the rate of energy 
loss through radiation. The rate of radiative cooling can be estimated using a steady state coronal 
equilibrium model by [Post et al. 19"77 , who provided cooling tables, Ca(T), for different elements. 



We note that radiation in lines is dominant in the total energy losses. In our experiments 
we vary the material composing the wires using either aluminium (Al), stainless steel (Fe) or 



tungsten (W) wires. Consideration of the Post et al. 1977 results show that for these materials 
Ca{T) has a single peak in the range T < lOOeV with the peak values 7 x 10^^^, 3 x 10^^^ and 
>> 5 X 10^^^ erg err? j s respectively for Al, Fe and W. Consideration of the cooling parameter 
X shows dc (X tc (X T/{niCA{T)). Thus we expect that x ^ 1 for strong shocks which form in the 
flow we produce in these experiments. Note that x decrease with increasing atomic number 
A. We expect the jets from tungsten wires, in particular, will be in the strongly radiative regime. 

In cylindrical wire array experiments the use of different materials, (and hence different 
cooling rates), resulted in a difference in equilibrium radius of the precursor plasma column. The 
equilibrium radius was significantly smaller for the high Z elements ( [Lebedev et al. 2001 , Chittenden 



et al. 2001). For conical array experiments the difference in cooling rates will lead to differences 



in jet collimation properties. In their description of the theory of converging conical flows. 



Canto et al. 1988| derived relations between the converging flow properties and the resulting 
jet characteristics. Of particular importance for the results presented here they found relations 
between the converging flow angle of incidence {9), jet radius (rj) and the half-opening angle a of 
the conical shock on the axis. Both rj and a depend on the inverse compression ratio {C, = Po/ Pi) 
behind the shock (where subscripts and 1 refer to upstream and downstream conditions). 



(l-C)-v/(l-C)2-4Ctan2^ 

tan a = ^ . 4-1 

itant^ 



tan a , . ^ 

tan d + tan a 
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In the latter expression yo is a measure of the width of the converging flow. The equations above 
demonstrate that as the compression ratio increases (and C, decreases) both the opening angle of 
the axial shock and the radius of the jet will decrease. Since post-shock compression will depend on 



the degree of radiative cooling, equations 4-1 and ^-2| predict that jets formed in our experiments 
will become more narrow as the atomic number associated with the plasma ions increases. Note 
the opening angle of the jet depends on the post-shock temperature and so we expect that it will 
also decrease with increased cooling (and increased atomic number). 

Typical images of jets formed in experiments using conical arrays of aluminium, stainless 
steel and tungsten wires are shown in Fig 3. These images demonstrate that the degree of jet 
collimation does strongly depend on the atomic number. Increasing coUimation is seen for the 
elements with higher A. The diameters of the jets at the end of the formation region, where the 
converging plasma flow is still present, decrease with increasing atomic number, being 3.5, 3 and 
2.5 mm for Al, Fe and W, respectively. For Al the jet is slightly diverging, while for stainless 
steel the observed opening angle is close to zero, and the jet radius remains constant over the 
propagation length of / ~ 15 mm. For tungsten, an apparent convergence of the jet occurs at the 
jet head. The radius of the jet decreases towards the tip and an opening angle becomes negative 
there. The shape of the tip may reflect the temporal history of the formation of conical standing 
shock however more study of this point is required. The observed dependence of jet collimation 
on atomic number is consistent with predictions of Canto's model and indicates that changes in 
jet radius/opening angle are due to a higher levels of radiative losses in high A plasmas. 

Experimental evidence of fast radiative cooling of the jet is also seen from gated soft x-ray 
images (not shown), filtered to transmit radiation in the interval 190 — 290 eV . The data show 
that emission from the jets is rapidly decreasing along the beam. The length of the emitting part 
of the jet is significantly smaller than that measured at the same time by laser probing for all 
materials tested in the experiments (Al, Fe and W). The characteristic length of the emission decay 
(from the base of the jet) is ~ 7 mm for Al, ?e ~ 5 mm for Fe and /e ~ 2 mm for W, while the 
characteristic lengths of the jets, seen in laser probing images at the same time, are Lj ~ 1.5 cm. 
For a 50 eV plasma we estimate that ~ 10% of the total radiation is in this passband. 

Finally we note that previous astrophysical studies of converging conical flows have 



been plagued by issues of stability. All studies to date, both numerical ( Tenorio-Tagle et 



al. 1988| , iMellema fc Frank 1997|) and analytical (|Cant6 et al. 1988| , [Frank, Bahck fc Livio 1996| ) 



have imposed cylindrical symmetry as a means to facilitate calculations. While these studies 
demonstrated that collimation would occur, it has not been clear if the collimation point (the 



conical shock) on the axis would be stable ( Garcia-Segura 1997 ). The experimental studies 



reported here shed light on this issue. It is important to recognize that the converging plasma 
flows formed in our experimental system have a substantial degree of modulation in azimuthal 
direction, determined by the number of wires in the array (16 in the present experiments). Thus 
one expects perturbations from axisymmetric converging conical flows with angular wavelength of 
5(f) ^ 23°. In spite of these modulations, the degree of the jet collimation seen in the experiments 
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was remarkably high. To further examine stabihty considerations we have performed additional 
experiments in which we deliberately introduced significant axial asymmetry in the fiow by 
removing 2 neighboring wires from the array. For the highly radiative case (W),a well-collimated 
jet was still generated though it emerged inclined from the axis of symmetry. These results 
suggests that a high degree of symmetry in conical flow models of proto-stellar and proto-planetary 
jets is not required if radiation losses are significant in the formation region. 

Finally we note that fiows formed by the present technique could be used to study propagation 



of astrophysical jets through ambient clouds, a situation observed in YSOs ( Reipurth et al. 1996 ). 
This opportunity is illustrated by the image shown in Fig. 4. In these proof-of-concept experiments, 
a thin plastic foil was installed perpendicular to the direction of the jet propagation. Soft x-ray 
emission from the conical standing shock on the array axis is absorbed in the foil converting its 
surface layer into plasma, which expands towards the jet. The gated soft x-ray image in Fig 4 
shows emission from the jet, which rapidly decays due to radiative cooling as the jet leaves the 
region of formation. X-ray emission is seen again from the region where the jet is decelerated by 
the foil as its bulk energy of motion is converted into thermal energy. This data is consistent with 
the conclusion that the thermal energy of the jet during propagation is significantly smaller than 
the directed kinetic energy, i.e., the Mach number of the jet is high. 



5. Conclusions and Astrophysical Implications 

We have presented results from the first of a series of experiments designed to create and 
study high mach number collimated plasma fiows. Our experiments were carried out on a pulsed 
power machine which drives IMA current through 16-wire array Z-pinch load. The wires are 
configured in a conical array such that as plasma is draw off the vaporized wires (by lorentz 
forces), it streams towards the axis at an oblique angle creating a converging conical fiow. As 
the fiow collides on the axis, a standing shock forms which redirects the plasma into a collimated 
beam (a jet.) 

Analysis of direct emission, laser probing and interferometric measurements of the jets reveal 
the following characteristics. 

• The jets are high mach number fiows with Vj ~ 200 km/ s and Mg > 20. 

• The jets are radiative in the sense that significant thermal energy losses to the emission of 
photons (predominantly in spectral lines). We estimate values of the cooling parameter to 
be X = dc/rj < 1 

• The fluid description of the plasma is valid {5 = X^fp/rj < 1), as is the inviscid assumption 
{Re >> 1) and the assumption that Heat conduction may be neglected to flrst order Pe > 1. 
We note that the values Pe are such that diffusive effects may produce smoothing of small 
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scale features. Based on the values of all the dimensionless numbers, a scaling between the 
behavior of the jets observed in the laboratory and astrophysical jets should be possible. 

The coUimation of the flows depends strongly in the efficacy of radiative cooling. As 
radiative cooling increases the jet radius and opening angle are observed to decrease. These 
measurements are consistent with the theory of Canto et al. 198^ . 



The formation of jets via converging conical flows is stable to azimuthal perturbations. 

Our experimental configuration allows direct experimental investigation of shock/cloud 
interactions. 



The results of these simulations are promising in that they may provide a means investigating 
fundamental issues relating to the collimation, propagation and stability of radiative plasma 
jets. The application of scaling relations between the laboratory experiments and astrophysical 
environments (particularly stellar jets i.e., YSOs and PNe) implies that these studies can be used 
to evaluate questions concerning the nature and evolution of collimated plasma flows in space. 
The experiments detailed in this study both complement and extend previous work in this area 
([Farely et al. 1999| , |Shi gemori et al. 2000| , ptone et al. 2000[ ) but provide an entirely different rout 



to creating radiative scalable plasma jets (Z-pinch vs. intense lasers). The larger size of the jets 
created in the Z-pinch systems, (scales of centimeters vs hundreds of microns), also facilitates 
imaging of the resultant plasma structures. 

The use of high energy density laboratory experiments to study cosmic plasma environments 



is, essentially, a new undertaking ( Remington et al. 1999 ). As such, the burden of proof will rest 



on the practitioners of the field to demonstrate the relevancy and effectiveness of this approach 
to astronomy. The critical issue is then to identify critical domains where the laboratory plasma 
studies can directly add information to our understanding of open astrophysical questions. The 
validation of astrophysical simulation codes using the experiments is one means of demonstrating 
relevance. It is hoped however that fundamentally new knowledge concerning astrophysical 
systems can be gleaned from these kinds of investigations. 

High mach number radiative jets has already shown some promise as being one domain 
where laboratory studies could prove relevant. In particular, fundamental issues relating to 
the nature of high mach number radiative jet/ambient interactions have never been explored 
experimentally. The mode structure of unstable jets also remains questions which would benefit 
from direct investigation. The studies described in this letter take an important first step 
towards demonstrating relevancy since they already appear to address one long standing issue in 
astrophysical jet studies: the stability of jet production via converging conical flows. While MHD 
effects may be the means by which winds are launched and collimated from stellar sources ( pudritzl 



1991 , Shu et al. 1994 , Blackman et al. 2001| ), numerical studies have shown that converging 



flows are likely to facilitate the creation of jets even when fields are present ( pardiner, Frank & 



Hartmann (2001)| ) 
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Thus while many issues remain to be studied 
as dynamically important magnetic fields, present 
relevant investigations of radiative plasma jets can 



in greater detail in future experiments, such 
results already suggest that astrophysically 
be carried out with Z-pinch laboratory studies. 
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Fig. 1. — Schematic of the experiment. 

Fig. 2. — Time-resolved laser probing images show well-coUimated plasma jet ejected from the 
tungsten wire array. Interferometry (left image) give electron density of the jet of 10^^ cm~^. Two 
schlieren images obtained in the same experiment with time delay of 30 ns give measurement of 
the jet velocity of 200 km/s. 

Fig. 3. — Laser probing images of plasma jets formed in Al, stainless steel and W wire arrays 
show that degree of coUimation increases for elements with higher atomic number, in which rate of 
radiative cooling is higher. 

Fig. 4. — Time-resolved soft x-ray image showing interaction of the jet with thin plastic foil. 
Emission from the jet rapidly decays due to radiative cooling as the jet leaves formation region and 
emission is seen again from the region where the jet is decelerated by the foil and kinetic energy is 
converted into thermal energy. 
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